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Abstract 

Using the technique of Pade Approximants for the correlation contributions of charged particles we are 
able to reproduce complex mathematical expressions by simple algebraic formulas. The Pade formulas are 
analytical expressions, which interpolate between certain density-temperature regions and are character- 
ized by exact asymptotics. We present a semi-relativistic description of the thermodynamics applicable 
for stellar interiors. Comparisons with equation-of-state data obtained by other calculational schemes 
are presented. 



1. Introduction 

The equation of state (EOS) formalism for astrophysical applications requires an accurate calcu- 
lation of the non-ideal effects in plasmas over a wide range of densities and temperatures. Many 
successful efforts have been made to provide EOS-tables [[[], [2], |3|, |5| ||, [?J which are applicable 
to describe thermodynamical properties of stellar matter. However, detailed comparisons of 
quantities composed of second-order derivatives of the free Helmholtz energy (e.g. specific heat, 
adiabatic temperature gradient) reveal considerable discontinuities (see, e.g., ||). We present 
an EOS-formalism which avoids such numerical deviations and is well suited for example to 
evaluate the thermodynamics of stellar interiors. 



2. Thermodynamical Functions 

We start with the thermodymamical description of a fully ionized plasma 
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where S = {F, G, PV, U, V/K^, $sV, CyT} and a = {/, g,p, u, efts, cy} symbolize the vari- 
ous thermodynamical functions, as the free Helmholtz energy F, free Gibbs energy G, pressure P, 
internal energy U, inverse isothermal compressibility 1/Kt, strain coefficient <3?S; and isochoric 
specific heat Cy. Quantities given by small letters are the corresponding potentials normalized 
to N e ^kBT. Terms labeled by id, x, and c denote contributions due to ideality, exchange, 
and correlation contributions of the plasma species a = e, i. The term labeled by cq takes into 
account quantum effects of ions. 

Furthermore we have to determine the isobaric specific heat Cp by the relation ^ 
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X p and xt are the density and temperature exponents in the EOS-formalism. The adiabatic 
temperature gradient defined by V ac j = (dlnT/dlnP)s (S is the entropy) is given by |8|, f?J 
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with the adiabatic exponents Ti = -y/PK^, T2 = 7/(7 — PKt^g) and T3 = 1 + 7c and the 
generalized Griineisen coefficient 7g = V&s/TCy. 

The ideal and exchange contributions, <rj, d and <j* e , are calculated in a relativistic approach, 
which covers the two approximations: a) arbitrary degeneracy (i/j) and weak relativity (A < 1) 
and b) strong degeneracy (ifi > 1) and arbitrary relativity (A = kT/mc 2 ). Details are given for 
free Helmholtz energy, free Gibbs energy, and pressure in ||. 

The correlation terms <7g e , erf;, and af e are taken into account by Pade approximants which 
are modified versions of those given in ||. The new Pade approximants are determined by a 
rearrangement of the quantum virial function which is dominant in the weakly coupled regime. 
Moreover the Pade formula for the ion-electron subsystem is extented by linear T-asymptotics 



for the classical region [1C]. The ionic quantum correction <r^ q are described by |11]. 

In this paper we present our new Pade approximant for the ion-ion correlation. Starting with 
the free energy the complete set of the quantities summarized by crS will be given. We mention 
that the Pade formulas for the quantities related by the first and second order derivatives 
of the free energy are constructed from the analytical limiting laws, which are known to be 
asymptotically exact. The intermediate region will be fitted by new parameters. Of course, 
the quality of this procedure with regard to the thermodynamic consistency is determined by 
comparisons with the numerical data produced from corresponding theories. The advantage 
is that the Pade approximants for our set of thermodynamic functions are expressed by their 
correct asymptotics, which cover a wide parameter-region. 

For the ionic subsystem we use the new formulas for the free Helmholtz- and Gibbs energy 
with Ti = (Z 5 / 3 )e 2 < 
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with 
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and the coefficients bo = l/\/3 ,61 = 3\/3/86q, 62 = 100, ^3 = bi/y/3, and to = 36o/2, t\ = 
5&i/3, t% = 62 1 t% = 463/3 . Eq. d|]) is based on the classical Ti < 1 result from Cohen and 
Murphy [Q| and for T; ^ 1 we take into account the Monte-Carlo fits for the free energy ey(ri) 
of the liquid and solid one-component plasma (OCP) from Stringfellow et al. fi"3|| . 

The internal energy, the compressibility, and the strain coefficient are determined by 
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Note, that a factor 1/3 is missing in the second term of 1/A^y given by Hansen 14]. 
For the isochoric specific heat we apply the Pade approximant 
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Fig. 1: Internal energy (left) and isochoric specific heat (right) of the OCP model. The 
short-dashed lines refer to the 1^ < 1 asymptotics from Cohen and Murphy [12] and the 
long-dashed line to the T; ^ 1 fit from Stringfellow et al. [JDJ. Left: The solid line refers 
to Eq. (0), the dots to Kahlbaum pq |, and the crosses to Chabrier and Potekhin []l^ ]. 
Right: The solid line refers to Eq. (||) and the dotted one to Hansen |H. 



with the coefficients qo = 3&o/4, q\ = 7b±, q% = 5&2> Q3 = 863. 
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3. Results and Discussion 



In Fig. | our new Pade approximants for the internal energy 
(|9|) are compared with other closed-form parametrizations p 
energy is already covered by the asymptotic theories given in 

detailed comparison between [|1| and the simplified fit formula from Chabrier and Potekhin [16 



and the isochoric specific heat 
[Tgl - Obviously, the internal 
13] over a wide Ti-region. A 




was carried out by Potekhin [17] recently. Comparisons between the internal energy derived from 
our earlier Pade approximants for the free energy || and a new fit formula from Chabrier 



and Potekhin [16] are shown in 16]. It is noteworthy, that our former Pade approximant for /j? 
proposed by Ebeling [18 contains ionic quantum corrections, which deliver in the classical limit 
the result of Abe [19]. The inclusion of quantum effects (see ||, the coefficient b\ in Eq. (65)) 
implies an additional temperature-dependence besides the pure T-dependence. Accordingly, the 
determination of thermal quantities such as internal energy, strain coefficient, or specific heat 
from the former Pade formula for f» must be carried out carefully. Non-thermal derivatives such 
as pressure did not suffer from this drawback as shown in ||. 

In our new calculational scheme all ionic quantum corrections are considered separately. As 
shown in Fig. ([!]) the classical interpolation formulas are well in agreement. 

Fig. ^illustrates the course of the adiabatic temperature gradient, V a d, calculated by our new 
Pade approximants and gives comparisons with numerical results from other authors |i~6| , 14, 20|. 



V a d is a highly sensitive quantity, because it depends on first and second-order derivatives of 
the model free Helmholtz energy as given by Eq. (Q). The low-density limes of V a< j = 0.25 (right 
panel) is a consequence of the photonic contribution in contrast to Fig. ^ (left panel). 

We have shown how by means of analytical expressions (see also ||, [H]] ) the thermodynamics 
of e.g. stellar interiors can be described over a wide range of densities and temperatures. 
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Fig. 2: Adiabatic temperature gradient V ac j vs. density for helium at T = 10 5 K 
without radiation (left) and for carbon at T = 10 7 K with radiation (right). Note, that 
the discontinuities at logp ~ 6 (left) and logp ~ 9 (right) are caused by the fluid-solid 
phase transition at T = 178 ]13|| . The solid curves refer to Eqs. (||), (||), (|l|) and the 
dashed-dotted to Q. Left: The dotted line refers to ||. For comparision is shown the 
ideality (long-dashed line) . Right The dashed line refers to |2lJ and the crosses to |^] . 



Our analytical description leads to a smooth course of V a d and avoids any numerical noise. 
Further comparisions are given in detail by p^ |. 
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